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Abstract 

The Allen Atlas of the adult mouse brain is a brain-wide, genome-wide data set that has been 
made available online, triggering a renaissance in neuroanatomy. In particular, it has been used to 
define brain regions in a computational, data-driven way, and to estimate the region-specificity of cell 
types characterized independently by their transcriptional activity. However, these results were based 
on one series of co-registered (coronal) ISH image series per gene, whereas the online ABA contains 
several image series per genes, including sagittal ones. Since the sagittal series cover mostly the left 
hemisphere, we can simulate the variability of results by repeatedly drawing a random image series 
for each gene and restricting the computation to the left hemisphere. This gives rise to an estimate 
of error bars on the results of computational neuroanatomy. This manuscript is a complement to 
the paper Computational neuroanatomy: mapping cell-type densities in the mouse brain, simulations 
from the Allen Brain Atlas prepared for the International Conference on Mathematical Modeling in 
Physical Sciences June 5-8, 2015, Mykonos Island, Greece. 
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Acronyms. ABA: Allen Brain Atlas; ARA: Allen Reference Atlas; ISH: in situ hybridization. 


Contents 


1 Introduction 

2 Methods [3] 

2.1 Gene expression data from the ABA. [3] 

2.2 Simulation of variability. [3] 

2.3 Correlation between the ABA and cell-type-specific microarray data. [4] 

2.4 Density estimates of cell types. [5] 

3 Results [6] 

3.1 Distributions of correlations give rise to peaks. [6] 

3.2 Peaks for the density of cell types are better separated than for correlations. [8] 

3.3 Ranking of cell types by stability of results . [10] 

4 Discussion an 

5 Tables: cell-type-specific transcriptomes: description, labeling and anatomical origin 1121 

6 Tables of of rankings of cell types by estimates of overlap with the coronal model 1151 

7 Cell-type-specific results 1161 


1 Introduction 

The Allen Brain Atlas (ABA, [T1 has renewed the gene-based approach to gene-expression studies 

in neuroscience by releasing voxelized, brain-wide ISH data for the entire genome of the mouse, which 
were co-registered to the Allen Reference Atlas (ARA, [24]). About 4,000 genes of special neurobiological 
interest were proritized. For these genes an entire brain was sliced coronally and processed (giving rise to 
the coronal ABA). For the rest of the genome the brain was sliced sagitally, and only the left hemisphere 
was processed (giving rise to the sagittal ABA). In the computational approaches of [191 130 El 123 IZ] ? 
only data from the coronal ABA were analyzed, in order to obtain brain-wide results. However, these 
brain-wide results revealed a large degree of left-right symmetry. Moreover, more than 90% of the genome 
was observed in mm to be expressed in the brain. There is therefore a strong need to incorporate data 
from the sagittal atlas into the analysis of the ABA. The website www.mouse-brain.org already caters 
to this need for each gene, because a query based on a gene name returns all the ISH image series for the 
corresponding gene (specifying whether it comes from sagittal or coronal sections), together with a bar 
diagram summarizing how the gene expression breaks up between the main regions of the brain. This 
enables a user of the atlas to compare the expression profiles coming from all the available data, on a 
gene-by-gene basis. The results are reproducible on a desktop computer using the Brain Gene Expression 
Analysis MATLAB toolbox [8]. 

The ABA has proven to be a powerful tool for systems biology, because co-registered gene-expression 
data can be studied collectively, thousands of genes at a time. Indeed the collective behaviour of gene- 
expression data is crucial for the analysis of m [2T1 123 5 which the brain-wide correlation between 
the ABA and cell-type-specific microarray data was studied. Strong heterogeneities of the correlation 
profiles were observed, allowing to guess, for instance that medium spiny neurons were extracted from 
the caudoputamen. Moreover, the region-specificity of cell types was estimated by linear regression with 
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positivity constraint. The model was fitted using the coronal atlas only, which allowed to obtain brain¬ 
wide results. However, this restriction implies the availability of only one ISH gene-expression profile per 
gene. This poses the problem of the error bars on the results of the model. 

In the present paper we shall extend the analysis to include ISH data from the sagittal atlas. These 
data are based on the sections of the left hemisphere only, which induces a restriction of neuroanatomical 
results to the left hemisphere, but allows to simulate the distribution of analysis results involving the 
collective behaviour of gene-expression profiles without restricting the number of genes. Since the ISH 
data in the sagittal atlas were co-registered to the voxelized ARA (at a resolution of 200 microns), they 
can be treated computationally in the same way as the coronal atlas. 

2 Methods 

2.1 Gene expression data from the ABA 

The digitized expression energies obtained by measuring the grey-scale intensities of ISH image series 
co-registered to the ARA mm can be arranged into matrices, in which rows correspond to voxels and 
columns to genes. At a spatial resolution of 200 microns, there are V = 49, 742 cubix voxels in the mouse 
brain). Consider a gene labeled g in the coronal atlas, and a voxel labeled v in the mouse brain. This 
gene comes with a number of gene-expression profiles in the ABA, call this number N{g). For most genes 
in the coronal atlas N{g) = 2, with one coronal profile and one sagittal profile, but there are cases with 
more data. Let us label these ISH profiles by integers between 1 and N{g) (with the label 1 corresponding 
to the coronal atlas, for definiteness). Let us denote by E^^\v^g) the gene-expression energy of gene g 
at voxel so that ^) is positive number corresponding to the gene-expression energy of gene g at 

voxel V in the coronal atla^ For higher values of the index /c, the entry E^^\v,g) is a positive number 
if voxel labeled v is covered in the tissue sections of sample labeled k, and is labeled as missing data 
otherwise (which is the case if the label k corresponds to a sagittal sample, and voxel v is in the right 
hemisphere). 

When working with all the image series in the atlas at the same time, we have to restrict to voxels 
that are in the intersection of the sets of voxels for which all genes have been processed A maximal- 
intensity projection of the sum of all gene-expression profile shown on Fig. which shows that the data 
are not homogeneous over the brain hemisphere on average, and that the left hemisphere indeed has a 
larger signal than the right hemisphere, because of the larger number of image series obtained in the ABA 
(the voxels with no data were treated as zero entries when performing the sum of gene expression energies). 


2.2 Simulation of variability 

For computational purposes we have the choice between N{g) expression profiles for gene labeled g. 
Hence, instead of just one voxel-by-gene matrix whose ^-th column contains , g)) , the ABA 

gives rise to a family of 11^=1 ^id) voxel-by-gene matrices, with voxels belonging to the left hemisphere. 
Any quantity computed from E^^^ can be recomputed from any of these matrices, thereby inducing a 
distribution for this quantity, from which we can estimate quantities such as mean, median, variance and 
density. This is a finite but prohibitively large number of computations (at least 2^, given that each gene 
in the coronal atlas has at least two image series) to be performed, so for practical purposes we have 
to take a Monte Carlo approach and to perform the computation for R random choices of data. The 
simulation of the distribution of a quantity Q depending on the atlas can be obtained by the following 
pseudo-code: 
for all i in [1..R] 

1. for all g in [1..G], choose an integer ni{g) in [l..N{g)] 

^This is the quantity denoted by E(v,g) in papers [To] [9] [8]|6] [20] 
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Figure 1: Maximal-intensity projections of the sum of all gene-expression profiles in the atlas (from both 
sagittal and coronal series). The left hemisphere has higher intensity then the right hemisphere due to 
the larger number of series covering it. However, the anatomical heterogeneity patterns (highlighting 
the hippocampus for instance), apart from intensity, seem to be left-right symmetric. However, there 
is a darker area on the left hemisphere that reveals missing sections from the sagittal series: indeed in 
this region the slices ae extremely small and fragile. The same issue is encountered in the most frontal 
sections of the olfactory buln in the coronal series. 


2. construct the matrix with entries E^i^{v^ g) = ^ g) 

3. compute the quantity Q using this matrix, call the result (3[^] 
end 

The larger R is, the more precise the estimates for the distribution of the quantity Q will be. 


2.3 Correlation between the ABA and cell-type-specific microarray data 

As an application, consider the correlation profile between the ABA and cell-type-specific transcriptome 
profiles collated in [25]. The microarray reads are arranged into a matrix C, whose rows correspond to cell 
types and whose columns correspond to genes (arranged in the same orders as in the matrix presentation 
of the ABA expression energies): C{t^g) is the microarray read for gene labeled g in the cell-type-specific 
sample labeled t, in the notations of [20|, were the index t take values between 1 and T = 64, and g 
takes values getween 1 and G = 2,131, which is the number of genes found in [20] to be found both in 
the coronal ABA and in all the microarray reads). If we represent the ABA by a matrix E^^^ of ISH data 
in the left hemisphere, obtained by drawing a random set in step 2 of the pseudo-code, the value of the 
brain-wide correlation at voxel v with cell type t reads: 

E^=i {C{t,g) - {E^i\{v.9) - ^YZ=iE\i]{v,g)) 

Corr[j](v,t) = ^ ^ (1) 

Y ELi {C{t, k)-^ fc)) Y!h=i (-E’w iv, h)-E 1 -E[i] {v, /i)) 

and V takes values corresponding to the left hemisphere. Having computed the family of correlation 
matrices (Corr[^])i<^<^, we can study the distribution of their t-th column for each t. For example we 
can compute the average correlation profile (where the average is performed over the index i labling the 
choice of image series): 

R 

(Corr(v, t)) = (2) 

i=l 

Moreover, we can estimate the dispersion of the correlation values for a cell type labeled t and a 
region labeled Vr in the ARA. For each random family of data (labeled by index i), we can compute the 
average correlation in the voxels belonging to region Vr according to the ARA (this time the average is 
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performed over voxels and the index of image series is not summed over): 

Corr^,[i](t) = -rfh XI Corr[i]{t>, i). 

I veVr 

For a fixed value of t, this computation furnishes us with one family of R numbers in the interval [—1,1] 
for each region in the ARA. For the correlation analyis of [20] to be stable under a change of animal and 
sectioning modality, these families of numbers should come from a probability density that presents a 
peak. 


2.4 Density estimates of cell types 

The linear model proposed in m estimates the density profiles of cell types characterized by the tran- 
scriptome profile, assuming the expression energy of each gene is proportional to the number of mRNA 
molecules at each voxel, and that the microarray read of each gene in a cell type is also proportional to 
the number of mRNAs for this gene in this cell type. The expression energy at voxel labeled v must be 
a a sum of cell-type-specific microarray reads, weighted by the density of each type at each voxel: 

-E(v, 5) = X Pt{v)C{t, g), (4) 

t 

where index t labels cell types. In [20|, we estimated the density profiles {pt{v)) for all voxels v in the 
mouse brain and for the cell types belonging to the data set collated in [25] . This estimation process was 
based on minimizing the difference between the l.h.s. and r.h.s. of Eq. [^ over all the possible positive 
coefficients p. This optimization procedure is deterministic, and relating the result to the mean density 
(decomposing the density into the sum of its mean and Gaussian noise) is a difficult problem in statistics 
(see [29]). 


Some error estimates on the value of pt{v) were obained in [20[ [22] using sub-sampling techniques, 
which involved mutilating the data by keeping only a random 10% of the coronal ABA, refitting the 
model, and repeating the operation. This induces a ranking of the cell types based on the stability of the 
results against sub-sampling. However, the set of genes on which the computation is based changes for 
each computation, and the fraction 10% is arbitrary (even though it is close to the fraction of the genome 
represented by our coronal data set). Having integrated the sagittal data into the atlas, we can now refit 
the model with the same sets of genes, only changing the set of image series from which they come, and 
we can do so in 11^=1 different ways. The only price we have to pay for this is the restriction of 
the results to the left hemisphere. This operation is just another application of the procedure outlined in 
the pseudo-code, with the role of the quantity Q played by the family of numbers (pt('^))i<t<T ’ 
voxels V in the left hemisphere. 

With the notations introduced above, we denote by (pt,[i]('^)) the density estimate obtained from the 
random matrix of ISH data : 


E[i] (v, 5) = X S') + Residual(v, g). 

t 


(5) 


We can group the voxels by region according to the ARA as we did for correlations in order to compare 
the results to classical neuroanatomy. The average density across random draws of image series for cell 
type lebeled t reads: 

(s>t(^)) = ^Xaw(^)- (6) 

i=l 

Since the number of cells of a given type in an extensive quantity, we compute the fraction of the total 
density contributed by each region, rather than the average density 0r, [i] (t) in region labeled r for sample 
labeled i and cell-type labeled t : 


4^r,[i] (^) 
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^^vGleft hemisphere 




X 

veVr 


(7) 
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(b) 

Figure 2: Heat maps of the average correlation between random image series from the ABA 
and cell-type-specific microarray data, (Corr(u, t)), defined in Eq. (a) For medium spiny neurons, 
t = 16. (b) For granule cells, t = 20. The sections were taken through the region in the ARA in which 
the correlation is highest on average (the boundary of this region defined by the ARA is outlined in green 
on the section). 


3 Results 

3.1 Distributions of correlations give rise to peaks 

After running simulations, we computed the average correlati8on profiles defined in Eq. This quantity 
is plotted is plotted as a heat map on Fig. for medium spiny neurons (t = 16) and granule cells 
(t = 20), from which it is clear that the average correlation profile presents a plateau inside the striatum 
for medium spiny neurons and inside the cerebellum for granule cells (which was also the case in [20] 
were the correlation profiles between the coronal atlas and cell types were analyzed). 

Choosing brain regions from the coarsest version of the ARA, we computed the regional averages of 
correlation profiles defined in Eq. The resulting families of numbers in the interval [—1,1] were 
subjected to density estimation (based on kernel methods in MATLAB). The estimated densities present 
peaks. Moreover, the peaks are localized at higher values of correlation for the regions that present 
striking correlation patterns. Results of density estimation for medium spiny neurons (t = 16) and 
granule cells (t = 20) are presented on Fig. In both cases the densities are peak-shaped, and the peak 
corresponding to the highest average correlation corresponds to the expected region, and is well decoupled 
from the other peaks: the peak corresponding to the striatum for the medium spiny neurons is centered 
at 0.29, and entirely supported in the interval [0.27,0.31], in which none of the quantities Corr^ [^](16) is 
found for values of r not corresponding to the striatum. These results induce reassuring bounds pointing 
at the stability of the results of [20l[2T], even though the dark areas in the coronal and axial projections 
of Fig. probably corresponds to an artifically low expression energy in the most lateral part of the left 
hemisphere, due to missing sections in sagittal series. This effect is probably bringing down the average 
correlation between the cortex and the medium spiny neurons, but it concerns only a small fraction of 
the cortex, and the correlation in the most lateral sections estimated from coronal series only is close to 
the cortical average. 
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Figure 3: Estimated probability densities of the correlation profiles agglomerated in a few regions of the 
coarsest versions of the ARA. All the region-specific densities give rise to peaks. The right-most peak is 
well-decoupled from the others in both cases, (a) For medium spiny neurons, t = 16; the peak centered 
at the highest value corresponds to the striatum, (b) For granule cells, t = 20; the peak centered at the 
highest value corresponds to the cerebellum. 
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(b) 

Figure 4: Heat maps of the average density of cell types in the left hemisphere, (Corr(u,t)), defined in 
Eq. 1^ (a) For medium spiny neurons, t = 16. (b) For granule cells, t = 20. They are visually very close 
to the heat maps of piQ and p 20 obtained in eq]. 

3.2 Peaks for the density of cell types are better separated than for correla¬ 
tions 

Having refitted the linear model (Eq. for each random draw of image series, we can plot the induced 
average density profiles for the same two cell types as above (see Fig. [^, which by eye gives a strong 
impression of the caudoputamen and of the granular layer of the cerebellum in the left hemisphere. More¬ 
over, we estimated the densities of the contribution of each region in the coarsest version of the ARA 
to the total density of each cell type (Eq. [^. The resulting densities are supported in the interval [0,1] 
by construction, and due to the fact the contributions from each region to the total density sum to 1 
in each sample, the righjt-most peaks have a tendency to be more clearly decoupled from the others 
than in the correlation analysis. In particular, we can read from Fig. that medium spiny neurons 
have 93(±3) percent supported in the striatum, and from Fig. similar numbers for granule cells in the 
cerebellum, without any region gathering more than 5 percent of the signal in any of the randoim samples. 

The two cases exposed in Figs. andfittingDistrs serve as proof of concept, because medium spiny 
neurons and granule cells are well-studied cell types, independently from the neurome approach of m , and 
they are expected to be strongly associated to the striatum and cerebellum respectively. However, some 
cell types may exhibit more complex neuroanatomical patterns of density, and breaking the estimated 
density according to the coarsest version of the ARA as in Eq. To study neuroanatomy in a purely 
data-driven way, we must not use any input from classical neuroanatomy. For a given cell type, we must 
compare the family of density profiles {p[i]^t)i<i<R fo the density profile pt estimated from the coronal 
atlas. A possible comparison, proposed in [20] to analyze the results of the sub-sampling simulations, 
involves the computation of the overlap between density profiles: 

Aht) ■- ^ > 0)/9[i],t(v), (8) 

which is the fraction of the total estimated density in the i-th draw supported by the coronal model. For 
each cell-type label t, the overlap T(.,t) is a random variable, whose distribution can be studied using 
the empirical cumulative distribution function (CDF) as follows: 

CDF{u,t) = ^ 10 6 [l..R],I{r,t) < w}|. (9) 
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Figure 5: Estimated probability densities of fractions of density agglomerated in a few regions of the 
coarsest versions of the ARA (see Eq. [^. All the region-specific densities give rise to peaks. The right¬ 
most peak is well-decoupled from the others in both cases, furthermore the other peaks are all centered 
close to zero (making several of them almost invisible), (a) Eor medium spiny neurons, t = 16; the peak 
centered at the highest value corresponds to the striatum, (b) Eor granule cells, t = 20; the peak centered 
at the highest value corresponds to the cerebellum. 
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Figure 6: A heat map of the CDFs of the random overlaps with the results of the model based on the 
coronal atlas. 

If we present the CDFs in matrix form, with one cell type per column, and rows corresponding to the 
value of the overlap, we can plot this matrix as a heat map, as on Fig. The more stable the model is 
for a cell type, the larger the dark area in the corresponding column is. This Figure has a much larger 
dark area than the analogous Figure in [20], in which the random step involved taking a random 10% of 
the genes and refitting the model. For instance, 10 cell types have an overlap of more than 95 % with 
the original model with probability 1, which was not the case for any of the cell types according to the 
sub-sampling procedure. The computational treatment of the full set of image series therefore reveals 
stronger stability properties of the linear model. 


3.3 Ranking of cell types by stability of results 

Let us recall the notations introduced in [22| to analyze the results of the simulated distribution of over¬ 
laps with the orioginal (coronal) model. 

Having simulated the distribution of the sub-sampled densities of all the T cell-type specific transcrip- 
tomes in our data set, we can estimate confidence thresholds in two ways, for a cell type labeled t. 

(1) Impose a threshold a in the interval [0,1] on the overlap with the density pt estimated in the 
linear model, and work out the probability pt,a of reaching that threshold from the sub-samples: 

Pt,a ■- > a) = i |s e [l..S'],X(s,t) > a|. (10) 

For a cell type labeled t, the distribution of the overlaps X(.,t) can be visualized using the cumulative 
distribution function CDF^ (in the space [0,1] of the values of the overlap between pt sub-sampled profiles 
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CDFt(a) = ||{sG[1..5],2:(s,i)<a}|. 

The value CDFt{u) is related to the probability defined in Eq. [^as follows: 

Pt,a = l-CDFt{a). ( 12 ) 

(2) Impose a threshold (3 in the interval [0,1] on the fraction of sub-samples, and work out which 
overlap Xthresh{t^ P) with the estimated density pt is reached by that fraction of the sub-samples. The 
threshold value of the intercept Tthresh{t^ P) is readily expressed in terms of the inverse of the cumulative 
distribution function: 

Ithresh{t,p) = CDYp{p). (13) 

The more stable the prediction pt is again sub-sampling, the more concentrated the values of X(.,t) 
are at high values (close to 1), the slower the take-off of the cumulative function CDFt is, the lower the 
value of CDFt((a) is, and the larger the probablity is (for a fixed value of a in [0,1]). 

For a fixed cell type labeled t, the values of p{t ,.) and Tthreshit^ •) can therefore be readily read off 
from a plot of the cumulative distribution function CDF (this plot is in the af3 plane in our notations). 
For the sake of visualization of results for all cell types, we constructed the matrix V, whose columns 
correspond to cell types, and whose rows correspond to values of the threshold a: 

;= CDFt(a), (14) 

where the index in the l.h.s. means that the cell types are ordered by decreasing order of overlap 

between the average sub-sampled profile and the predicted profile. If the entries of the matrix V are 
plotted as a heat map (see Fig. [^, the hot colors will be more concentrated in the left-most part of the 
image. For a fixed column, the more concentrated the hot colors are in the heat map, the more stable 
the corresponding cell type is. 

4 Discussion 

The present work allows to bring the power of the online facilities of the ABA to the desktop for a col¬ 
lective analysis of gene-expression energies. The example of the analysis of correlation with microarray 
data sets shows the power of the ABA as a tool of analysis for other data sets. This goes beyond the 
correlation structure of the ABA itself, and its relation to the classical neuroanatomy of the ARA, that 
can elready be investigated online using the Anatomic Gene Expression Atlas (AGFA, see [2]). 

When the coronal ABA was used to analyze other data sets, in particular microarray data as in [20] . 
it was assumed that it reflects an average collective expression of the genome in the adult mouse brain, 
showing little sensitivity to dynamics and to animal-to-animal variations, even though each expression 
profile comes from a different animal. The fact that correlation profiles allowed to guess the anatomical 
origin of many cell types in [20] provided a self-consistent check of this assumption, but incorporating 
the different ISH image series into the computational analysis allows to control the numerical errors in 
the correlation profiles. 

Moreover, the Monte Garlo approach we took to simulate the variability of region-specificity is much 
more consistent than the sub-sempling approach taken in [20j. The numerical results show a greater 
stability than in the sub-sampling approach on average. 

It should be emphasized that the computations in the present paper only serve as examples, bringing 
more accurate answers to the most current questions posed by the analysis of the ABA. Ultimately, 
the distribution of any quantity based on the atlas can be computed exactly and studied at each voxel, 
without the need to agglomerate voxels into neuroanatomical regions. 
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5 Tables: cell-type-specific transcriptomes: description, label¬ 
ing and anatomical origin 

For each of the cell-type-specific samples analyzed in this note, the following two tables give a brief 
description of the cell type, the region from which the samples were extracted according to the coarsest 
version of the Allen Reference Atlas, and the finest region to which it can be assigned according to the 
data provided in the studies [niiiiiiiniiiziiiiiiiniiiaiii]. The indices in the first columns of the tables 
are the ones refered to as t. 


12 


Index 

Description 

Region in the ARA 
(’bigl2’) 

Finest label in the ARA 

1 

Purkinje Cells 

Gerebellum 

Gerebellar cortex 

2 

Pyramidal Neurons 

Gerebral cortex 

Primary motor area; Layer 5 

3 

Pyramidal Neurons 

Gerebral cortex 

Primary somatosensory area; Layer 5 

4 

A9 Dopaminergic Neurons 

Midbrain 

Substantia nigra_ compact part 

5 

AlO Dopaminergic Neurons 

Midbrain 

Ventral tegmental area 

6 

Pyramidal Neurons 

Gerebral cortex 

Gerebral cortex; Layer 5 

7 

Pyramidal Neurons 

Gerebral cortex 

Gerebral cortex; Layer 5 

8 

Pyramidal Neurons 

Gerebral cortex 

Gerebral cortex; Layer 6 

9 

Mixed Neurons 

Gerebral cortex 

Gerebral cortex 

10 

Motor Neurons, Midbrain Cholinergic Neurons 

Midbrain 

Peduncolopontine nucleus 

11 

Cholinergic Projection Neurons 

Pallidum 

Palhdum_ ventral region 

12 

Motor Neurons, Cholinergic Interneurons 

Medulla 

Spinal cord 

13 

Cholinergic Neurons 

Striatum 

Striatum 

14 

Interneurons 

Gerebral cortex 

Gerebral cortex 

15 

Drdl+ Medium Spiny Neurons 

Striatum 

Striatum 

16 

Drd2+ Medium Spiny Neurons 

Striatum 

Striatum 

17 

Golgi Cells 

Gerebellum 

Gerebellar cortex 

18 

Unipolar Brush cells (some Bergman 
Glia) 

Gerebellum 

Gerebellar cortex 

19 

Stellate Basket Cells 

Gerebellum 

Gerebellar cortex 

20 

Granule Cells 

Gerebellum 

Gerebellar cortex 

21 

Mature Oligodendrocytes 

Gerebellum 

Gerebellar cortex 

22 

Mature Oligodendrocytes 

Gerebral cortex 

Gerebral cortex 

23 

Mixed Oligodendrocytes 

Gerebellum 

Gerebellar cortex 

24 

Mixed Oligodendrocytes 

Gerebral cortex 

Gerebral cortex 

25 

Purkinje Cells 

Gerebellum 

Gerebellar cortex 

26 

Neurons 

Gerebral cortex 

Gerebral cortex 

27 

Bergman Glia 

Gerebellum 

Gerebellar cortex 

28 

Astroglia 

Gerebellum 

Gerebellar cortex 

29 

Astroglia 

Gerebral cortex 

Gerebral cortex 

30 

Astrocytes 

Gerebral cortex 

Gerebral cortex 

31 

Astrocytes 

Gerebral cortex 

Gerebral cortex 

32 

Astrocytes 

Gerebral cortex 

Gerebral cortex 

33 

Mixed Neurons 

Gerebral cortex 

Gerebral cortex 

34 

Mixed Neurons 

Gerebral cortex 

Gerebral cortex 

35 

Mature Oligodendrocytes 

Gerebral cortex 

Gerebral cortex 

36 

Oligodendrocytes 

Gerebral cortex 

Gerebral cortex 

37 

Oligodendrocyte Precursors 

Gerebral cortex 

Gerebral cortex 


Table 1: Anatomical origin of the cell-type-specific samples (I). 
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Index 

Description 

Region in the ARA 
(’bigl2’) 

Finest label in the ARA 

38 

Pyramidal Neurons, Callosally projecting, P3 

Gerebral cortex 

Gerebral cortex 

39 

Pyramidal Neurons, Callosally projecting, P6 

Gerebral cortex 

Gerebral cortex 

40 

Pyramidal Neurons, Callosally projecting, P14 

Gerebral cortex 

Gerebral cortex 

41 

Pyramidal Neurons, Corticospinal, P3 

Gerebral cortex 

Gerebral cortex 

42 

Pyramidal Neurons, Corticospinal, P6 

Gerebral cortex 

Gerebral cortex 

43 

Pyramidal Neurons, Corticospinal, P14 

Gerebral cortex 

Gerebral cortex 

44 

Pyramidal Neurons, Corticotectal, P14 

Gerebral cortex 

Gerebral cortex 

45 

Pyramidal Neurons 

Gerebral cortex 

Gerebral cortex. Layer 5 

46 

Pyramidal Neurons 

Gerebral cortex 

Gerebral cortex. Layer 5 

47 

Pyramidal Neurons 

Gerebral cortex 

Primary somatosensory area; Layer 5 

48 

Pyramidal Neurons 

Gerebral cortex 

Prelimbic area and Infralimbic area; Layer 

5 (Amygdala) 

49 

Pyramidal Neurons 

Hippocampal 

region 

Ammon’s Horn; Layer 6B 

50 

Pyramidal Neurons 

Gerebral cortex 

Primary motor area 

51 

Tyrosine Hydroxylase Expressing 

Pons 

Pontine central gray 

52 

Purkinje Cells 

Gerebellum 

Gerebellar cortex 

53 

Glutamatergic Neuron (not well defined) 

Gerebral cortex 

Gerebral cortex; Layer 6B (Amygdala) 

54 

GABAergic Interneurons, VIP+ 

Gerebral cortex 

Prelimbic area and Infralimbic area 

55 

GABAergic Interneurons, VIP+ 

Gerebral cortex 

Primary somatosensory area 

56 

GABAergic Interneurons, SST+ 

Gerebral cortex 

Prelimbic area and Infralimbic area 

57 

GABAergic Interneurons, SST-|- 

Hippocampal 

region 

Ammon’s Horn 

58 

GABAergic Interneurons, PV-f- 

Gerebral cortex 

Prelimbic area and Infralimbic area 

59 

GABAergic Interneurons, PV-t- 

Thalamus 

Dorsal part of the lateral geniculate com¬ 
plex 

60 

GABAergic Interneurons, PV+, P7 

Gerebral cortex 

Primary somatosensory area 

61 

GABAergic Interneurons, PV+, PIO 

Gerebral cortex 

Primary somatosensory area 

62 

GABAergic Interneurons, PV+, P13-P15 

Gerebral cortex 

Primary somatosensory area 

63 

GABAergic Interneurons, PV+, P25 

Gerebral cortex 

Primary somatosensory area 

64 

GABAergic Interneurons, PV+ 

Gerebral cortex 

Primary motor area 


Table 2: Anatomical origin of the cell-type-specific samples (II). 
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Cell type 

Granule Cells 
Purkinje Cells 

Motor Neurons, Cholinergic Interneurons 
Pyramidal Neurons 
Pyramidal Neurons 

Glutamatergic Neuron (not well defined) 

Mature Oligodendrocytes 

Drd2+ Medium Spiny Neurons 

Tyrosine Hydroxylase Expressing 

Pyramidal Neurons, Callosally projecting, P14 

Pyramidal Neurons 

Astroglia 

Pyramidal Neurons, Corticotectal, P14 
Pyramidal Neurons 
GABAergic Interneurons, PV+ 

GABAergic Interneurons, SST + 


Purkinje Cells 

GABAergic Interneurons, SST + 

GABAergic Interneurons, PV+ 

Mature Oligodendrocytes 
Pyramidal Neurons 
AlO Dopaminergic Neurons 
Cholinergic Neurons 
Pyramidal Neurons 
GABAergic Interneurons, VIP + 

Motor Neurons, Midbrain Cholinergic Neurons 
A9 Dopaminergic Neurons 
Oligodendrocyte Precursors 
Mature Oligodendrocytes 

Pyramidal Neurons, Callosally projecting, P3 


Index t 

X{t), (%) 

Pt,0.759 (%) 

'^thn 

20 

99.5 

100 

99.2 

1 

99.4 

100 

99.2 

12 

99.2 

100 

99.2 

47 

99.1 

100 

99.2 

49 

99.1 

100 

99.2 

53 

99.1 

100 

99.6 

35 

99 

100 

98.8 

16 

98.9 

100 

98.8 

51 

98.9 

100 

98.8 

40 

98.3 

100 

98.4 

46 

98 

100 

98.4 

28 

97.1 

100 

98.4 

44 

96.9 

100 

97.7 

6 

92.4 

99.6 

95.3 

64 

92.2 

100 

93 

57 

91.7 

100 

93.4 

31 

91.1 

99.4 

93.4 

52 

89.9 

93.9 

95.3 

56 

88.2 

95.7 

93 

59 

87 

85.8 

94.5 

21 

87 

100 

88.3 

48 

86.2 

98.9 

88.7 

5 

85.7 

84.1 

93.4 

13 

84.7 

96.9 

87.5 

45 

82.8 

97.7 

85.2 

55 

81.6 

75 

89.5 

10 

80.7 

94.1 

82.8 

4 

80.1 

74.8 

87.5 

37 

79.6 

67.1 

91.8 

22 

79 

78.4 

81.6 

38 

78.5 

70.9 

95.7 

32 

78.3 

65.5 

93 


Table 3: Cell types ranked by average overlap with the coronal model (I). 


Tables of of rankings of cell types by estimates of overlap with 
the coronal model 




Cell type 

Stellate Basket Cells 
GABAergic Interneurons, PV+, P7 
Mixed Neurons 

Pyramidal Neurons, Corticospinal, P14 
Pyramidal Neurons 
GABAergic Interneurons, VIP + 

Astrocytes 

Pyramidal Neurons, Callosally projecting, P6 

Mixed Neurons 

Oligodendrocytes 

Pyramidal Neurons 

Drdl+ Medium Spiny Neurons 

Golgi Cells 

Pyramidal Neurons, Corticospinal, P6 
GABAergic Interneurons, PV+, P25 
GABAergic Interneurons, PV+ 

Unipolar Brush cells (some Bergman Glia) 
GABAergic Interneurons, PV+, PIO 
Astroglia 

Pyramidal Neurons 
Pyramidal Neurons 
Neurons 
Mixed Neurons 
Bergman Glia 
Pyramidal Neurons 

GABAergic Interneurons, PV+, P13-P15 

Cholinergic Projection Neurons 

Interneurons 

Mixed Oligodendrocytes 

Mixed Oligodendrocytes 

Purkinje Cells 

Pyramidal Neurons, Corticospinal, P3 


Index t 

j(t), (%) 

Pt,0.759 (%) 

'^thn 

19 

77 

64.5 

80.1 

60 

76.3 

57.1 

89.5 

34 

72.8 

50.3 

83.2 

43 

69.9 

48.2 

87.1 

7 

69 

18.5 

73.4 

54 

67.4 

15 

72.7 

30 

63.5 

9.9 

69.1 

39 

56.2 

40.5 

89.1 

33 

56.2 

21.5 

71.9 

36 

54.3 

8.5 

63.7 

50 

44.4 

9.7 

62.9 

15 

42.8 

-0.4 

47.7 

17 

40.4 

36.8 

99.2 

42 

32 

9.3 

48.8 

63 

30.4 

-0.3 

39.1 

58 

26.2 

0.6 

35.5 

18 

25 

-0.4 

31.3 

61 

17.9 

3.1 

29.3 

29 

14.1 

-0.4 

18 

8 

12.1 

0.1 

14.1 

3 

9.6 

8 

0.4 

26 

4 

0.1 

5.9 

9 

3.8 

3.1 

0.4 

27 

1.1 

-0.4 

0.4 

2 

0.3 

-0.4 

0.4 

62 

0 

-0.4 

0.4 

11 

0 

22.2 

71.9 

14 

0 

22.2 

71.9 

23 

0 

22.2 

71.9 

24 

0 

22.2 

71.9 

25 

0 

22.2 

71.9 

41 

0 

22.2 

71.9 


Table 4: Cell types ranked by average overlap with the coronal model (II). 

7 Cell-type-specific results 
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Figure 7: Heat map of the mean correlation profile for cell type labeled 1. 
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Figure 8: Distribution profiles of the correlations between transcriptome profile of cell type labeled 1 
across random sets of ISH expression energy profiles. 



Figure 9: Heat map of the mean density profile for cell type labeled 1. 
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Figure 10: Distribution profiles of the density profile of cell type labeled 1 across random sets of ISH 
expression energy profiles. 
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Figure 11: Heat map of the mean correlation profile for cell type labeled 2. 
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Figure 12: Distribution profiles of the correlations between transcriptome profile of cell type labeled 2 
across random sets of ISH expression energy profiles. 



Figure 13: Heat map of the mean density profile for cell type labeled 2. 
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Figure 14: Distribution profiles of the density profile of cell type labeled 2 across random sets of ISH 
expression energy profiles. 
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Figure 15: Heat map of the mean correlation profile for cell type labeled 3. 
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Figure 16: Distribution profiles of the correlations between transcriptome profile of cell type labeled 3 
across random sets of ISH expression energy profiles. 



Figure 17: Heat map of the mean density profile for cell type labeled 3. 
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Figure 18: Distribution profiles of the density profile of cell type labeled 3 across random sets of ISH 
expression energy profiles. 
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Figure 19: Heat map of the mean correlation profile for cell type labeled 4. 
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Figure 20: Distribution profiles of the correlations between transcriptome profile of cell type labeled 4 
across random sets of ISH expression energy profiles. 
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Figure 21: Heat map of the mean density profile for cell type labeled 4. 
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Figure 22: Distribution profiles of the density profile of cell type labeled 4 across random sets of ISH 
expression energy profiles. 
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Figure 23: Heat map of the mean correlation profile for cell type labeled 5. 
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Figure 24: Distribution profiles of the correlations between transcriptome profile of cell type labeled 5 
across random sets of ISH expression energy profiles. 



Figure 25: Heat map of the mean density profile for cell type labeled 5. 



■ Basic cell groups and regions 

- Cerebral cortex 
Olfactory areas 

■Hippocampal region 

■ Retrohippocampal region 
Striatum 

■Pallidum 

Thalamus 

- Hypothalamus 
Midbrain 

■Pons 

■Medulla 

Cerebellum 


0 0.5 

Fraction of total density 


Figure 26: Distribution profiles of the density profile of cell type labeled 5 across random sets of ISH 
expression energy profiles. 
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Figure 27: Heat map of the mean correlation profile for cell type labeled 6. 
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Figure 28: Distribution profiles of the correlations between transcriptome profile of cell type labeled 6 
across random sets of ISH expression energy profiles. 



Figure 29: Heat map of the mean density profile for cell type labeled 6. 
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Figure 30: Distribution profiles of the density profile of cell type labeled 6 across random sets of ISH 
expression energy profiles. 
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Figure 31: Heat map of the mean correlation profile for cell type labeled 7. 
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Figure 32: Distribution profiles of the correlations between transcriptome profile of cell type labeled 7 
across random sets of ISH expression energy profiles. 



Figure 33: Heat map of the mean density profile for cell type labeled 7. 



Basic cell groups and regions 

- Cerebral cortex 

Olfactory areas 
^—Hippocampal region 

Retrohippocampal region 

Striatum 

Pallidum 

Thalamus 

- Hypothalamus 

Midbrain 

Pons 

^—Medulla 

Cerebellum 


Fraction of total density 


Figure 34: Distribution profiles of the density profile of cell type labeled 7 across random sets of ISH 
expression energy profiles. 
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Figure 35: Heat map of the mean correlation profile for cell type labeled 8. 
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Figure 36: Distribution profiles of the correlations between transcriptome profile of cell type labeled 8 
across random sets of ISH expression energy profiles. 



Figure 37: Heat map of the mean density profile for cell type labeled 8. 


to 

c 

0) 

Q 



Basic cell groups and regions 

- Cerebral cortex 

Olfactory areas 
^—Hippocampal region 

Retrohippocampal region 

Striatum 

Pallidum 

Thalamus 

- Hypothalamus 

Midbrain 

Pons 

^—Medulla 

Cerebellum 


Fraction of total density 


Figure 38: Distribution profiles of the density profile of cell type labeled 8 across random sets of ISH 
expression energy profiles. 
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Figure 39: Heat map of the mean correlation profile for cell type labeled 9. 
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Figure 40: Distribution profiles of the correlations between transcriptome profile of cell type labeled 9 
across random sets of ISH expression energy profiles. 



Figure 41: Heat map of the mean density profile for cell type labeled 9. 
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Figure 42: Distribution profiles of the density profile of cell type labeled 9 across random sets of ISH 
expression energy profiles. 
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Figure 43: Heat map of the mean correlation profile for cell type labeled 10. 
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Figure 44: Distribution profiles of the correlations between transcriptome profile of cell type labeled 10 
across random sets of ISH expression energy profiles. 




Figure 45: Heat map of the mean density profile for cell type labeled 10. 


0.25 



■ Basic cell groups and regions 

- Cerebral cortex 
Olfactory areas 

■Hippocampal region 

■ Retrohippocampal region 
Striatum 

■Pallidum 

Thalamus 

- Hypothalamus 
Midbrain 

■Pons 

■Medulla 

Cerebellum 


0.2 0.4 0.6 0.8 

Fraction of total density 


Figure 46: Distribution profiles of the density profile of cell type labeled 10 across random sets of ISH 
expression energy profiles. 
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Figure 47: Heat map of the mean correlation profile for cell type labeled 11. 
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Figure 48: Distribution profiles of the correlations between transcriptome profile of cell type labeled 11 
across random sets of ISH expression energy profiles. 



Figure 49: Heat map of the mean density profile for cell type labeled 11. 
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Figure 50: Distribution profiles of the density profile of cell type labeled 11 across random sets of ISH 
expression energy profiles. 
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Figure 51: Heat map of the mean correlation profile for cell type labeled 12. 
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Figure 52: Distribution profiles of the correlations between transcriptome profile of cell type labeled 12 
across random sets of ISH expression energy profiles. 



Figure 53: Heat map of the mean density profile for cell type labeled 12. 
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Figure 54: Distribution profiles of the density profile of cell type labeled 12 across random sets of ISH 
expression energy profiles. 
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Figure 55: Heat map of the mean correlation profile for cell type labeled 13. 
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Figure 56: Distribution profiles of the correlations between transcriptome profile of cell type labeled 13 
across random sets of ISH expression energy profiles. 
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Figure 57: Heat map of the mean density profile for cell type labeled 13. 
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Figure 58: Distribution profiles of the density profile of cell type labeled 13 across random sets of ISH 
expression energy profiles. 
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Figure 59: Heat map of the mean correlation profile for cell type labeled 14. 
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Figure 60: Distribution profiles of the correlations between transcriptome profile of cell type labeled 14 
across random sets of ISH expression energy profiles. 



Figure 61: Heat map of the mean density profile for cell type labeled 14. 
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Figure 62: Distribution profiles of the density profile of cell type labeled 14 across random sets of ISH 
expression energy profiles. 
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Figure 63: Heat map of the mean correlation profile for cell type labeled 15. 
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Figure 64: Distribution profiles of the correlations between transcriptome profile of cell type labeled 15 
across random sets of ISH expression energy profiles. 
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Figure 65: Heat map of the mean density profile for cell type labeled 15. 
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Figure 66: Distribution profiles of the density profile of cell type labeled 15 across random sets of ISH 
expression energy profiles. 
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Figure 67: Heat map of the mean correlation profile for cell type labeled 16. 



Figure 68: Distribution profiles of the correlations between transcriptome profile of cell type labeled 16 
across random sets of ISH expression energy profiles. 



Figure 69: Heat map of the mean density profile for cell type labeled 16. 



Figure 70: Distribution profiles of the density profile of cell type labeled 16 across random sets of ISH 
expression energy profiles. 
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Figure 71: Heat map of the mean correlation profile for cell type labeled 17. 
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Figure 72: Distribution profiles of the correlations between transcriptome profile of cell type labeled 17 
across random sets of ISH expression energy profiles. 



Figure 73: Heat map of the mean density profile for cell type labeled 17. 
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Figure 74: Distribution profiles of the density profile of cell type labeled 17 across random sets of ISH 
expression energy profiles. 
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Figure 75: Heat map of the mean correlation profile for cell type labeled 18. 
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Figure 76: Distribution profiles of the correlations between transcriptome profile of cell type labeled 18 
across random sets of ISH expression energy profiles. 



Figure 77: Heat map of the mean density profile for cell type labeled 18. 
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Figure 78: Distribution profiles of the density profile of cell type labeled 18 across random sets of ISH 
expression energy profiles. 
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Figure 80: Distribution profiles of the correlations between transcriptome profile of cell type labeled 19 
across random sets of ISH expression energy profiles. 
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Figure 81: Heat map of the mean density profile for cell type labeled 19. 
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Figure 82: Distribution profiles of the density profile of cell type labeled 19 across random sets of ISH 
expression energy profiles. 
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Figure 83: Heat map of the mean correlation profile for cell type labeled 20. 
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Figure 84: Distribution profiles of the correlations between transcriptome profile of cell type labeled 20 
across random sets of ISH expression energy profiles. 



Figure 85: Heat map of the mean density profile for cell type labeled 20. 



Figure 86: Distribution profiles of the density profile of cell type labeled 20 across random sets of ISH 
expression energy profiles. 
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Figure 87: Heat map of the mean correlation profile for cell type labeled 21. 
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Figure 88: Distribution profiles of the correlations between transcriptome profile of cell type labeled 21 
across random sets of ISH expression energy profiles. 



Figure 89: Heat map of the mean density profile for cell type labeled 21. 
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Figure 90: Distribution profiles of the density profile of cell type labeled 21 across random sets of ISH 
expression energy profiles. 
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Figure 91: Heat map of the mean correlation profile for cell type labeled 22. 
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Figure 92: Distribution profiles of the correlations between transcriptome profile of cell type labeled 22 
across random sets of ISH expression energy profiles. 



Figure 93: Heat map of the mean density profile for cell type labeled 22. 
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Figure 94: Distribution profiles of the density profile of cell type labeled 22 across random sets of ISH 
expression energy profiles. 
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Figure 95: Heat map of the mean correlation profile for cell type labeled 23. 
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Figure 96: Distribution profiles of the correlations between transcriptome profile of cell type labeled 23 
across random sets of ISH expression energy profiles. 



Figure 97: Heat map of the mean density profile for cell type labeled 23. 
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Figure 98: Distribution profiles of the density profile of cell type labeled 23 across random sets of ISH 
expression energy profiles. 
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Figure 99: Heat map of the mean correlation profile for cell type labeled 24. 
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Figure 100: Distribution profiles of the correlations between transcriptome profile of cell type labeled 24 
across random sets of ISH expression energy profiles. 



Figure 101: Heat map of the mean density profile for cell type labeled 24. 


w 

c 

<u 

Q 


1 |-^-!- 1 -!- 

0.5 -.....- 

0 : i : : 

-0.5 -.....- 

^^^^- 

0 0.2 0.4 0.6 0.8 1 


Fraction of total density 


Basic cell groups and regions 

- Cerebral cortex 

Olfactory areas 
^—Hippocampal region 

Retrohippocampal region 

Striatum 

Pallidum 

Thalamus 

- Hypothalamus 

Midbrain 

Pons 

^—Medulla 

Cerebellum 


Figure 102: Distribution profiles of the density profile of cell type labeled 24 across random sets of ISH 
expression energy profiles. 
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Figure 103: Heat map of the mean correlation profile for cell type labeled 25. 
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Figure 104: Distribution profiles of the correlations between transcriptome profile of cell type labeled 25 
across random sets of ISH expression energy profiles. 



Figure 105: Heat map of the mean density profile for cell type labeled 25. 
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Figure 106: Distribution profiles of the density profile of cell type labeled 25 across random sets of ISH 
expression energy profiles. 
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Figure 107: Heat map of the mean correlation profile for cell type labeled 26. 


c 

(U 

O 



Basic cell groups and regions 

- Cerebral cortex 

Olfactory areas 
^—Hippocampal region 

Retrohippocampal region 
— Striatum 

- Pallidum 

Thalamus 

- Hypothalamus 

Midbrain 

- Pons 

- Medulla 

Cerebellum 


Figure 108: Distribution profiles of the correlations between transcriptome profile of cell type labeled 26 
across random sets of ISH expression energy profiles. 



Figure 109: Heat map of the mean density profile for cell type labeled 26. 
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Figure 110: Distribution profiles of the density profile of cell type labeled 26 across random sets of ISH 
expression energy profiles. 
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Figure 111: Heat map of the mean correlation profile for cell type labeled 27. 
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Figure 112: Distribution profiles of the correlations between transcriptome profile of cell type labeled 27 
across random sets of ISH expression energy profiles. 



Figure 113: Heat map of the mean density profile for cell type labeled 27. 
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Figure 114: Distribution profiles of the density profile of cell type labeled 27 across random sets of ISH 
expression energy profiles. 
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Figure 115: Heat map of the mean correlation profile for cell type labeled 28. 
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Figure 116: Distribution profiles of the correlations between transcriptome profile of cell type labeled 28 
across random sets of ISH expression energy profiles. 



Figure 117: Heat map of the mean density profile for cell type labeled 28. 
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Figure 118: Distribution profiles of the density profile of cell type labeled 28 across random sets of ISH 
expression energy profiles. 
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Figure 119: Heat map of the mean correlation profile for cell type labeled 29. 
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Figure 120: Distribution profiles of the correlations between transcriptome profile of cell type labeled 29 
across random sets of ISH expression energy profiles. 



Figure 121: Heat map of the mean density profile for cell type labeled 29. 
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Figure 122: Distribution profiles of the density profile of cell type labeled 29 across random sets of ISH 
expression energy profiles. 
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Figure 123: Heat map of the mean correlation profile for cell type labeled 30. 
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Figure 124: Distribution profiles of the correlations between transcriptome profile of cell type labeled 30 
across random sets of ISH expression energy profiles. 



Figure 125: Heat map of the mean density profile for cell type labeled 30. 
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Figure 126: Distribution profiles of the density profile of cell type labeled 30 across random sets of ISH 
expression energy profiles. 
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Figure 127: Heat map of the mean correlation profile for cell type labeled 31. 
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Figure 128: Distribution profiles of the correlations between transcriptome profile of cell type labeled 31 
across random sets of ISH expression energy profiles. 



Figure 129: Heat map of the mean density profile for cell type labeled 31. 
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Figure 130: Distribution profiles of the density profile of cell type labeled 31 across random sets of ISH 
expression energy profiles. 
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Figure 131: Heat map of the mean correlation profile for cell type labeled 32. 
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Figure 132: Distribution profiles of the correlations between transcriptome profile of cell type labeled 32 
across random sets of ISH expression energy profiles. 



Figure 133: Heat map of the mean density profile for cell type labeled 32. 
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Figure 134: Distribution profiles of the density profile of cell type labeled 32 across random sets of ISH 
expression energy profiles. 
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Figure 135: Heat map of the mean correlation profile for cell type labeled 33. 
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Figure 136: Distribution profiles of the correlations between transcriptome profile of cell type labeled 33 
across random sets of ISH expression energy profiles. 



Figure 137: Heat map of the mean density profile for cell type labeled 33. 
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Figure 138: Distribution profiles of the density profile of cell type labeled 33 across random sets of ISH 
expression energy profiles. 
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Figure 139: Heat map of the mean correlation profile for cell type labeled 34. 
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Figure 140: Distribution profiles of the correlations between transcriptome profile of cell type labeled 34 
across random sets of ISH expression energy profiles. 



Figure 141: Heat map of the mean density profile for cell type labeled 34. 
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Figure 142: Distribution profiles of the density profile of cell type labeled 34 across random sets of ISH 
expression energy profiles. 
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Figure 143: Heat map of the mean correlation profile for cell type labeled 35. 
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Figure 144: Distribution profiles of the correlations between transcriptome profile of cell type labeled 35 
across random sets of ISH expression energy profiles. 



Figure 145: Heat map of the mean density profile for cell type labeled 35. 
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Figure 146: Distribution profiles of the density profile of cell type labeled 35 across random sets of ISH 
expression energy profiles. 
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Figure 147: Heat map of the mean correlation profile for cell type labeled 36. 
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Figure 148: Distribution profiles of the correlations between transcriptome profile of cell type labeled 36 
across random sets of ISH expression energy profiles. 



Figure 149: Heat map of the mean density profile for cell type labeled 36. 
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Figure 150: Distribution profiles of the density profile of cell type labeled 36 across random sets of ISH 
expression energy profiles. 
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Figure 151: Heat map of the mean correlation profile for cell type labeled 37. 
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Figure 152: Distribution profiles of the correlations between transcriptome profile of cell type labeled 37 
across random sets of ISH expression energy profiles. 
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Figure 153: Heat map of the mean density profile for cell type labeled 37. 
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Figure 154: Distribution profiles of the density profile of cell type labeled 37 across random sets of ISH 
expression energy profiles. 
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Figure 155: Heat map of the mean correlation profile for cell type labeled 38. 
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Figure 156: Distribution profiles of the correlations between transcriptome profile of cell type labeled 38 
across random sets of ISH expression energy profiles. 



Figure 157: Heat map of the mean density profile for cell type labeled 38. 
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Figure 158: Distribution profiles of the density profile of cell type labeled 38 across random sets of ISH 
expression energy profiles. 
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Figure 159: Heat map of the mean correlation profile for cell type labeled 39. 
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Figure 160: Distribution profiles of the correlations between transcriptome profile of cell type labeled 39 
across random sets of ISH expression energy profiles. 



Figure 161: Heat map of the mean density profile for cell type labeled 39. 
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Figure 162: Distribution profiles of the density profile of cell type labeled 39 across random sets of ISH 
expression energy profiles. 
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Figure 163: Heat map of the mean correlation profile for cell type labeled 40. 
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Figure 164: Distribution profiles of the correlations between transcriptome profile of cell type labeled 40 
across random sets of ISH expression energy profiles. 



Figure 165: Heat map of the mean density profile for cell type labeled 40. 
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Figure 166: Distribution profiles of the density profile of cell type labeled 40 across random sets of ISH 
expression energy profiles. 


56 






































Figure 167: Heat map of the mean correlation profile for cell type labeled 41. 
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Figure 168: Distribution profiles of the correlations between transcriptome profile of cell type labeled 41 
across random sets of ISH expression energy profiles. 



Figure 169: Heat map of the mean density profile for cell type labeled 41. 
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Figure 170: Distribution profiles of the density profile of cell type labeled 41 across random sets of ISH 
expression energy profiles. 
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Figure 171: Heat map of the mean correlation profile for cell type labeled 42. 
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Figure 172: Distribution profiles of the correlations between transcriptome profile of cell type labeled 42 
across random sets of ISH expression energy profiles. 



Figure 173: Heat map of the mean density profile for cell type labeled 42. 
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Figure 174: Distribution profiles of the density profile of cell type labeled 42 across random sets of ISH 
expression energy profiles. 
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Figure 175: Heat map of the mean correlation profile for cell type labeled 43. 
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Figure 176: Distribution profiles of the correlations between transcriptome profile of cell type labeled 43 
across random sets of ISH expression energy profiles. 



Figure 177: Heat map of the mean density profile for cell type labeled 43. 
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Figure 178: Distribution profiles of the density profile of cell type labeled 43 across random sets of ISH 
expression energy profiles. 
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Figure 179: Heat map of the mean correlation profile for cell type labeled 44. 
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Figure 180: Distribution profiles of the correlations between transcriptome profile of cell type labeled 44 
across random sets of ISH expression energy profiles. 



Figure 181: Heat map of the mean density profile for cell type labeled 44. 
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Figure 182: Distribution profiles of the density profile of cell type labeled 44 across random sets of ISH 
expression energy profiles. 


60 













































Sagittal proj. Coronal proj. Axial proj. Plane of section Section 



Figure 183: Heat map of the mean correlation profile for cell type labeled 45. 
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Figure 184: Distribution profiles of the correlations between transcriptome profile of cell type labeled 45 
across random sets of ISH expression energy profiles. 
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Figure 185: Heat map of the mean density profile for cell type labeled 45. 
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Figure 186: Distribution profiles of the density profile of cell type labeled 45 across random sets of ISH 
expression energy profiles. 
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Figure 187: Heat map of the mean correlation profile for cell type labeled 46. 
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Figure 188: Distribution profiles of the correlations between transcriptome profile of cell type labeled 46 
across random sets of ISH expression energy profiles. 



Figure 189: Heat map of the mean density profile for cell type labeled 46. 
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Figure 190: Distribution profiles of the density profile of cell type labeled 46 across random sets of ISH 
expression energy profiles. 
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Figure 191: Heat map of the mean correlation profile for cell type labeled 47. 
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Figure 192: Distribution profiles of the correlations between transcriptome profile of cell type labeled 47 
across random sets of ISH expression energy profiles. 




Figure 193: Heat map of the mean density profile for cell type labeled 47. 
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Figure 194: Distribution profiles of the density profile of cell type labeled 47 across random sets of ISH 
expression energy profiles. 
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Figure 195: Heat map of the mean correlation profile for cell type labeled 48. 
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Figure 196: Distribution profiles of the correlations between transcriptome profile of cell type labeled 48 
across random sets of ISH expression energy profiles. 



Figure 197: Heat map of the mean density profile for cell type labeled 48. 



■ Basic cell groups and regions 

- Cerebral cortex 
Olfactory areas 

■Hippocampal region 

■ Retrohippocampal region 
Striatum 

■Pallidum 

Thalamus 

- Hypothalamus 
Midbrain 

■Pons 

■Medulla 

Cerebellum 


0.1 0.2 0.3 0.4 0.5 

Fraction of total density 


Figure 198: Distribution profiles of the density profile of cell type labeled 48 across random sets of ISH 
expression energy profiles. 
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Figure 199: Heat map of the mean correlation profile for cell type labeled 49. 
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Figure 200: Distribution profiles of the correlations between transcriptome profile of cell type labeled 49 
across random sets of ISH expression energy profiles. 



Figure 201: Heat map of the mean density profile for cell type labeled 49. 
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Figure 202: Distribution profiles of the density profile of cell type labeled 49 across random sets of ISH 
expression energy profiles. 
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Figure 203: Heat map of the mean correlation profile for cell type labeled 50. 


c 

(U 

Q 


0.04 
0.035 
0.03 
0.025 
0.02 
0.015 
0.01 
0.005 
0 

-0.2 -0.15 -0.1 -0.05 0 0.05 0.1 0.15 0.2 

Correlation 



Basic cell groups and regions 

- Cerebral cortex 

Olfactory areas 
^—Hippocampal region 

Retrohippocampal region 
— Striatum 

- Pallidum 

Thalamus 

- Hypothalamus 

Midbrain 

- Pons 

- Medulla 

Cerebellum 


Figure 204: Distribution profiles of the correlations between transcriptome profile of cell type labeled 50 
across random sets of ISH expression energy profiles. 



Figure 205: Heat map of the mean density profile for cell type labeled 50. 
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Figure 206: Distribution profiles of the density profile of cell type labeled 50 across random sets of ISH 
expression energy profiles. 
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Figure 207: Heat map of the mean correlation profile for cell type labeled 51. 
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Figure 208: Distribution profiles of the correlations between transcriptome profile of cell type labeled 51 
across random sets of ISH expression energy profiles. 
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Figure 209: Heat map of the mean density profile for cell type labeled 51. 
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Figure 210: Distribution profiles of the density profile of cell type labeled 51 across random sets of ISH 
expression energy profiles. 
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Figure 211: Heat map of the mean correlation profile for cell type labeled 52. 
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Figure 212: Distribution profiles of the correlations between transcriptome profile of cell type labeled 52 
across random sets of ISH expression energy profiles. 



Figure 213: Heat map of the mean density profile for cell type labeled 52. 
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Figure 214: Distribution profiles of the density profile of cell type labeled 52 across random sets of ISH 
expression energy profiles. 
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Figure 215: Heat map of the mean correlation profile for cell type labeled 53. 
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Figure 216: Distribution profiles of the correlations between transcriptome profile of cell type labeled 53 
across random sets of ISH expression energy profiles. 



Figure 217: Heat map of the mean density profile for cell type labeled 53. 
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Figure 218: Distribution profiles of the density profile of cell type labeled 53 across random sets of ISH 
expression energy profiles. 
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Figure 219: Heat map of the mean correlation profile for cell type labeled 54. 
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Figure 220: Distribution profiles of the correlations between transcriptome profile of cell type labeled 54 
across random sets of ISH expression energy profiles. 
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Figure 221: Heat map of the mean density profile for cell type labeled 54. 
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Figure 222: Distribution profiles of the density profile of cell type labeled 54 across random sets of ISH 
expression energy profiles. 
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Figure 223: Heat map of the mean correlation profile for cell type labeled 55. 
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Figure 224: Distribution profiles of the correlations between transcriptome profile of cell type labeled 55 
across random sets of ISH expression energy profiles. 



Figure 225: Heat map of the mean density profile for cell type labeled 55. 
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Figure 226: Distribution profiles of the density profile of cell type labeled 55 across random sets of ISH 
expression energy profiles. 
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Figure 227: Heat map of the mean correlation profile for cell type labeled 56. 
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Figure 228: Distribution profiles of the correlations between transcriptome profile of cell type labeled 56 
across random sets of ISH expression energy profiles. 



Figure 229: Heat map of the mean density profile for cell type labeled 56. 
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Figure 230: Distribution profiles of the density profile of cell type labeled 56 across random sets of ISH 
expression energy profiles. 
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Figure 231: Heat map of the mean correlation profile for cell type labeled 57. 
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Figure 232: Distribution profiles of the correlations between transcriptome profile of cell type labeled 57 
across random sets of ISH expression energy profiles. 




Figure 233: Heat map of the mean density profile for cell type labeled 57. 
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Figure 234: Distribution profiles of the density profile of cell type labeled 57 across random sets of ISH 
expression energy profiles. 
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Figure 235: Heat map of the mean correlation profile for cell type labeled 58. 
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Figure 236: Distribution profiles of the correlations between transcriptome profile of cell type labeled 58 
across random sets of ISH expression energy profiles. 



Figure 237: Heat map of the mean density profile for cell type labeled 58. 


w 

c 

<u 

Q 


1 |-^-!- 1 -!- 

0.5 -.....- 

0 : i : : 

- 0.5 -.....- 

^^^^- 

0 0.2 0.4 0.6 0.8 1 


Fraction of total density 


Basic cell groups and regions 

- Cerebral cortex 

Olfactory areas 
^—Hippocampal region 

Retrohippocampal region 

Striatum 

Pallidum 

Thalamus 

- Hypothalamus 

Midbrain 

Pons 

^—Medulla 

Cerebellum 


Figure 238: Distribution profiles of the density profile of cell type labeled 58 across random sets of ISH 
expression energy profiles. 
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Figure 239: Heat map of the mean correlation profile for cell type labeled 59. 
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Figure 240: Distribution profiles of the correlations between transcriptome profile of cell type labeled 59 
across random sets of ISH expression energy profiles. 




Figure 241: Heat map of the mean density profile for cell type labeled 59. 
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Figure 242: Distribution profiles of the density profile of cell type labeled 59 across random sets of ISH 
expression energy profiles. 
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Figure 243: Heat map of the mean correlation profile for cell type labeled 60. 
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Figure 244: Distribution profiles of the correlations between transcriptome profile of cell type labeled 60 
across random sets of ISH expression energy profiles. 



Figure 245: Heat map of the mean density profile for cell type labeled 60. 
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Figure 246: Distribution profiles of the density profile of cell type labeled 60 across random sets of ISH 
expression energy profiles. 
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Figure 247: Heat map of the mean correlation profile for cell type labeled 61. 
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Figure 248: Distribution profiles of the correlations between transcriptome profile of cell type labeled 61 
across random sets of ISH expression energy profiles. 



Figure 249: Heat map of the mean density profile for cell type labeled 61. 
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Figure 250: Distribution profiles of the density profile of cell type labeled 61 across random sets of ISH 
expression energy profiles. 
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Figure 251: Heat map of the mean correlation profile for cell type labeled 62. 


c 

(U 

O 



- 0.15 


- 0.05 


0.05 


Basic cell groups and regions 

- Cerebral cortex 

Olfactory areas 
^—Hippocampal region 

Retrohippocampal region 
— Striatum 

- Pallidum 

Thalamus 

- Hypothalamus 

Midbrain 

- Pons 

- Medulla 

Cerebellum 


Correlation 


Figure 252: Distribution profiles of the correlations between transcriptome profile of cell type labeled 62 
across random sets of ISH expression energy profiles. 



Figure 253: Heat map of the mean density profile for cell type labeled 62. 
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Figure 254: Distribution profiles of the density profile of cell type labeled 62 across random sets of ISH 
expression energy profiles. 
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Figure 255: Heat map of the mean correlation profile for cell type labeled 63. 
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Figure 256: Distribution profiles of the correlations between transcriptome profile of cell type labeled 63 
across random sets of ISH expression energy profiles. 



Figure 257: Heat map of the mean density profile for cell type labeled 63. 
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Figure 258: Distribution profiles of the density profile of cell type labeled 63 across random sets of ISH 
expression energy profiles. 
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Figure 259: Heat map of the mean correlation profile for cell type labeled 64. 
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Figure 260: Distribution profiles of the correlations between transcriptome profile of cell type labeled 64 
across random sets of ISH expression energy profiles. 



Figure 261: Heat map of the mean density profile for cell type labeled 64. 
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Figure 262: Distribution profiles of the density profile of cell type labeled 64 across random sets of ISH 
expression energy profiles. 
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